Scanning Tunneling Microscope Operating as a Spin-diode 
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We theoretically investigate spin-polarized transport in a system composed of a ferromagnetic 
Scanning Tunneling Microscope (STM) tip coupled to an adsorbed atom (adatom) on a host surface. 
Electrons can tunnel directly from the tip to the surface or via the adatom. Since the tip is 
ferromagnetic and the host surface (metal or semiconductor) is non-magnetic we obtain a spin-diode 
effect when the adatom is in the regime of single occupancy. This effect leads to an unpolarized 
current for direct bias {V > 0) and polarized current for reverse {V < 0) bias voltages, if the 
tip is nearby the adatom. Within the nonequilibrium Keldysh technique we analyze the interplay 
between the lateral displacement of the tip and the intra adatom Coulomb interaction on the spin- 
diode effect. As the tip moves away from the adatom the spin-diode effect vanishes and the currents 
become polarized for both V > and V < 0. We also find an imbalance between the up and down 
spin populations in the adatom, which can be tuned by the tip position and the bias. Finally, due 
to the presence of the adsorbate on the surface, we observe spin-resolved Friedel oscillations in the 
current, which reflects the oscillations in the calculated LDOS of the subsystem surface-l- adatom. 



I. INTRODUCTION 



The Scanning Tunneling Microscope (STM) has al- 
lowed huge advances in condensed matter physics. On 
one hand it serves as a powerful tool to manipulate mat- 
ter on a single atomic scale; ^'^ on the other it is used as 
a probe to the topology of metallic and semiconductor 
surfaces.^ An impressive early example of such control is 
the quantum corral, assembled by moving atom-by-atom 
on a metallic surface."''^ 

In the fascinating field of spintronics,^ STM was re- 
cently used to manipulate individual Co atoms adsorbed 
on a template of Mn.^ It was possible, for instance, to 
determine the spin direction of the individual Co atoms. 
STM was also applied to study the interactions between 
isolated Mn acceptors and the influence of the surface 
on the impurity properties in diluted magnetic semi- 
conductors, e.g. Mn-dopcd GaAs.®'^ More interesting, 
spin-polarized STM, sensitive to surface magnetization,^*^ 
has been used to map the morphology and the den- 
sity of states of single magnetic structures'^ and mag- 
netic quantum dots.''^ STM has also been employed in 
the investigation and identification of promising molec- 
ular switches, which could be used in future nanoscale 
circuits. "'^'^"'^ In the context of quantum information,'^ 
STM was used to measure electron spin relaxation times 
of individual atoms adsorbed on a surface with nanosec- 
ond time resolution. '^''^ More recently a new type of 
Scanning Probe Microscope was demonstrated using ul- 
tracold atoms. Interestingly, the conventional solid tip 
is replaced by a gas of ultracold rubidium atoms, which 
increases the spatial resolution of the microscope. All 
these applications highlight the importance of STM to 
the development of nano- engineered systems for spintron- 
ics and spin-based quantum information processing. 

As for strong correlated systems, STM has revealed 
new facets of the Kondo physics, such as the Fano-Kondo 
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FIG. 1. Ferromagnetic Scanning Tunneling Microscope (FM- 
STM) tip coupled to a host surface with an adatom. The 
matrix elements t\2, ii3 and t23 represent the couplings tip- 
adatom, tip-surface and adatom-surface, respectively. The 
tip-adatom lateral distance is denoted by R. 



resonance on the differential conductance when a sin- 
gle magnetic adatom (e.g., Co) is placed on a metallic 
surface. If the tip is ferromagnetic the Fano-Kondo 
line shape becomes spin-dependent'^' and the setup can 
be used as a powerful spin filter. 

Here we study spin-dependent transport in a system 
composed of a ferromagnetic (FM) STM tip coupled to 
both an adsorbed atom and a host non-magnetic (NM) 
surface. This geometry resembles a junction NM-QD- 
FM, where "QD" stands for quantum dot. In the present 
system, however, we have one additional ingredient: the 
tip can move laterally away from the adatom. It is well 
known that the NM-QD-FM system gives rise to unpo- 
larized current for direct bias {V > 0) and polarized cur- 
rent for reverse bias {V < 0) when the dot is singly oc- 
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cupicd (Coulomb blockade forbids double oecupancy).^^ 
This rectification of the current polarization is the so 
called spin-diode effect. In this work we investigate the 
interplay between the Coulomb blockade and the lateral 
displacement of the tip on the spin-diode effect. 

Our calculation is based on the Keldysh nonequilib- 
rium technique. By performing a self-consistent calcula- 
tion we determine the adatom occupation and its magne- 
tization as a function of the tip lateral position R, Fig. 1. 
We find that the adatom becomes magnetized when the 
tip approaches it; this magnetization switches sign as the 
bias is reversed. Moreover, we calculate the spin-resolved 
current in both single and double occupancy regimes of 
the adatom and find that it can depend strongly on the 
tip position. In particular, the spin-diode effect emerges 
when the tip is closer to the adatom and the charging en- 
ergy is large enough to allow for only a single electron in 
the adatom. As the tip moves away from the adatom 
the currents become equally polarized for both biases 
thus resulting in a suppression of the spin-diode behav- 
ior. More interesting, we observe spin-resolved Friedel 
oscillations^^"^^ in the current as the lateral distance tip- 
adatom R increases. These are due to the presence of 
the adatom on the surface, and reflect the oscillations 
in the local density of states (LDOS) of the subsystem 
surface-l- adatom. 

The paper is organized as follows. In Sec. II we present 
a detailed description of the formulation used to compute 
the spin-resolved currents and the spin populations. We 
divide this section into three subsections: in A we deter- 
mine an expression for the spin-resolved currents and the 
local density of states of the subsystem adatom-surface 
when the tip is far away from the adatom. In B, the 
calculation for the current is performed in the presence 
of the tip, and in C we summarize the numerical tech- 
nique and the parameters used in the simulations. We 
present our results and discussions in Sec. III. Section 
IV summarizes the main ideas of our work. 



II. FORMULATION 

The system we study is composed of a FM-STM tip 
over an adsorbate on a host surface. Fig. 1. Electrons 
can tunnel from the tip to the substrate or vice- versa via 
two possible ways: (i) direct tunneling tip-to-surface or 
(ii) tunneling via the localized state of the adatom. The 
system Hamiltonian is 

H = Hi + H2 + Hs + Hi2 + Hr3 + H23, (1) 

where Hi corresponds to the tip (i = I), the adatom 
(i = 2) or the surface (i = 3), and reads 

Hi^^ ek,<jc[^o-Cki<T + Si2Uh^ni, (2) 

kiCr 

where for i = 1(3), is the wave number for electrons 
in the tip (host) and the label cr stands for the electron 



spin components t and J,. Here ek^o- is the energy of the 
state kicr and Cki^ (c^ ,^) annihilates (creates) an elec- 
tron in the quantum state kicr. We consider a Stoner- 
like ferromagnetic dispersion ekia = h^kf/2m + crA for 
the tip, with m being the free electron mass and A the 
usual Stoner parameter, ^^'^^ and a free electron disper- 
sion Eka = f?k^/2m for the surface. For the adatom, 
1 = 2, we consider only a single spin-degenerate energy 
level, EkaCT = fcr- In this case the index k2 simply denotes 
the adatom level. The second term in H2 accounts for 
the Coulomb interaction U in the adatom. 

The coupling terms in Eq. (1) can be written as 

Hij= 51 (*y^k.^^k,'T +ty4^„Ck.a), (3) 

kikj (T 

where is the coupling parameter between subsystems 
i and j; ^12, ii3 and i23 account for the tunnelings 
tip-adatom, tip-surface and adatom-surface, respectively. 
When a bias voltage is applied these transfer terms drive 
the system out of equilibrium. 

Next we consider a real space formulation for the spin- 
resolved current. This is particularly convenient since we 
are interested in looking at Friedel oscillations on the sur- 
face. As we shall see later on, this formulation is equiv- 
alent to a formulation in the k space. 

A. Non-resonant transport 

For simplicity let us first consider the transport regime 
in which the direct coupling between the tip and the 
adatom is negligible (non- resonant transport), which is 
valid for large enough tip-adatom lateral distances. The 
Hamiltonian of the system in this case reduces to 

H = Hi+H2 + H3 + Hi3 + H23. (4) 

The electrical current for spin a between the tip and 
the surface can be calculated from the definition'^^ 

^-e{N^)^-te{[H,Nn), (5) 

with e the electron charge (e > 0) and A^f the total 
number operator given by 

Af = J dnvl/f (ri,i)*?(ri,t), (6) 

where 4'J'(ri,t) and \1/J'*(ri,f) are quantum field opera- 
tors for the electrons in the tip. In Eq. (5) and through- 
out the paper we assume ft = 1 . 

The quantity {0{t)) defines the nonequilibrium aver- 
age value of a physical observable denoted by the opera- 
tor 0{t), and it is given by'^^''^^ 

(0(i))=Tr[pO(t)], (7) 

where p is the thermal equilibrium density matrix, p = 
{TTe~^^")~^e^^^" , with Hq being the Hamiltonian con- 
taining only the Hi terms in Eq. (4), and 0{t) is in the 
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Hcisenbcrg picture, i.e., its time-dependence is governed 
by the full Hamiltonian of Eq. (4). 

The only non- vanishing term in the commutator of Eq. 
(5) is [i?i3, N^]. For electrons with spin cr, the tip-surface 
coupling can be written as^'^ 



Hiz — 



\j J rfridr3[T(ri,r3)*f (ri,<)*^(r3,i)+M, 

(8) 

where ^(ri, r^) is a matrix element that accounts for the 
coupling between the tip and the surface, and vl/|^(r3,t) 
is the quantum field operator for electrons in the surface. 
Calculating [His, N^] and using the result in Eq. (5) we 
find for the spin-resolved current 



I^ = ie J J dridr3[T(ri,r3)(vI'J (n, t)vl/-(r3, i)) - 

r*(ri,r3)(*f (r3,i)*?(ri,t))].(9) 
Defining the lesser Green function 

G<(r3, is; ri, ii) = i{^f (ri, ti)*^ (rg, h)), (10) 
we can rewrite the current as 



/f = 2eRc|y y dridr3T(ri,r3)G<(r3,i;ri,0|. (11) 

We now aim at determining G^(r3, t; ri,t) in Eq. (11). 
To this end, we use the noncquilibrium Kcldysh formal- 
ism. Similarly to the equilibrium case, here we introduce 
an ordered Green function 

G,(r3,T3;ri,Ti) ==-i(rc*^(r3,T3)*f (ri,Ti)), (12) 

with the r's defined, however, on a contour G in the com- 
plex plane. The operator Tc, called contour-ordering op- 
erator, orders the operators according to the position of 
their time arguments on the contour. From the contour- 
ordered Green function we can obtain the lesser G^, 
greater G^, retarded G^ and advanced G° Green func- 
tions, which are directly linked to the observables. 

To obtain G^, G^, GJ^ and G°, we first determine the 
equation of motion for the Green function in Eq. (12), 

a \ f 
^dn ~ 2rn) '^'^('^3, ^3; ri, n) = - / (ir^jT(r^j, ri) x 

G<,(r3,T3;rJj,T;^), 



(13) 



or in the integral form 



G<^(r3,r3;ri,Ti) = J J dr[dr'^ J dfG^(r3, T3; 4, f) x 

r(r;„r;)g,(r;,f;ri,ri), (14) 



where g^, is the free-electron Green function of the tip and 
the time integral is over the contour G ; then we perform 
an appropriate analytical continuation. This procedure 
consists essentially in replacing the contour integral over 



T in Eq. (13) by a real time integral over Here we follow 
the Langreth procedure.^"* For the lesser Green function 
G5 we have 



G<(r3,t3;ri,ti) 



dr[dr'^ I dt x 



G<(r3,t3;r^,i)T(r^,rl),9^(r;,t;ri,<i)] . (15) 

In the above equation and g'^ correspond to the 
analytically continued free-electron advanced and lesser 
Green functions of the tip, respectively. Throughout the 
paper wc use lower case to denote the free-electron Green 
functions of the tip, the adatom and the surface. We 
note that G^ (rs, ^3; ri, ti) is coupled to GJ^(r3, ^3; rg, f) 
and also to G^(r3, ^3; rg, t). To completely determine 
G^(r3, ^3; ri, ti) we then need to perform an iterative 
process and obtain a system of equations for the Green 
functions G^ and G^. 

Substituting Eq. (15) into Eq. (11) we obtain 




/i" = 2eRe<; / / / / j dvidv:idv'^dv'^dtT{vi,V:i)T{v'^,v[) x 
[Gl{vs,t-v'^,i)g<{v',,i-vi,t) + 

G<(r3,t;r;5,%^(r'i,f;ri,t)]|. (16) 

Performing a Fourier transform in the time coordinate 
we find 



" / J J J dvidv^dv'.dv'^ X 



drirfrsdr'^drg x (17) 
r(ri , r3 )r(r^5 , ) [G^ (r3 , , c.)g< (r^ , n , c.) + 
G<(r3,r^„c.)5^(r;,ri,c.)]|. 

Additionally, Fourier transforming g"^ and g'^ results in 

I" 

J 2tt 

e'-^i-^^T{vi,v:,)e''''<T{v'^,v[) (18) 
[G;(r3, r^j, c.)g< ^(c) + G<(r3, v'^,u:)gl^^{u)] | . 

Now we assume that the coupling between the tip and 
the surface is locaP'^''^^ (point source), i.e., T(ri,r3) = 
T°(ri)i5(r3 — R), r°(ri) being a proportionality factor. 
Integrating over the surface coordinates we find 

I- = 2e j ^Rci^ j j dridr'ie-'''i '-irO(ri) x 

I ki 

e^'^-^iTO*(r;)[G;(R,R,a.)g<,(c.) + 

G<(R,R,w)g£^,(c.)]l, (19) 
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and then over ri and r[ we arrive at 
/f=2e|^Re|^|tVPx 



(20) 



[G;(R,R,w)g<,M + G<(R,R,c.).g£,.H] 



where f?3kj is the Fourier transform of T"(ri). It is con- 
venient to perform a Fourier transform on GJJ.(R, R, w) 
and G^(R, R, w). This results in 



I'^ = 2e 



/ ^R4EE*i3k,k3W;,x (21) 

lk3,ki ki 



density of states per spin, P3(w) = J^k^^i^ ~ ^k,,) = 
m/mrh^ ~ 1/2D = pg, with n the 2D electronic density 
and D the band half width. 

In this case and in the zero temperature limit, the cur- 
rent reads 



(25) 



where V is the bias voltage. 

As we shall see in the next section, in the presence of 
the adatom, similarly to the case of an impurity on the 
surface of a metal, the LDOS shows Friedel oscillations; 
these affect the current that oscillates around the value 
given by Eq. (25). 



[GLk;,..H.<.M + G< ^,.(^)5£,.(^)]|, 

where tiakikg = ^iski^*'^'^^-'^^ what follows we as- 
sume t^ski ^ constant, i.e., independent of ki. 
As we have mentioned before, Eq. (21) could have 
been directly derived from the k space Hamiltonian 
Hij, = Eki,k3*i3Cki<TCk3<T + tJaCkscrCkia, with ti3 = 
tj'ae*''^'^. Observe that this is completely equivalent to 
the real space Hamiltonian Eq. (8), with a point source, 
T(ri,r3)=rO(ri)<5(r3-R). 



1. Spatially resolved transmission coefficient 

Equation (20) can alternatively be written as 

Ii=e y'^[2^|4lW(^)]*{G<(R,R,a.)+(22) 

/i(c.)[G;(R,R,c.)-G^(R,R,a;)]}, 

where /i(w) and pia{<^) = X^ki '^('^ ^ ^kicr) are, respec- 
tively, the Fermi function and the density of states of 
the tip. If the tip is weakly coupled to the surface, we 
can calculate the current in the lowest order of |ti3p. 
This means that the Green functions G^(R, R, w) and 
G^(R, R, cj) are assumed tip-decoupled, i.e., here they 
are local Green functions for the subsystem adatom- 
surface only. In other words, the tip behaves as a probe, 
not affecting the local density of states. This allows us 
to write the current as 



2^ 



[/i(a.)-/3(c^)]T,(R,L.), (23) 



where f^ioj) is the Fermi function of the surface and the 
transmission coefficient is defined by 

r,(R,c.) = [27r|<?3|2pi,(,.)] [-2Im{G;(R,R,c.)}] 

= [2n\t%\''pi,{Lo)\ [2^p,(R,c^)], (24) 

with p<^(R,a;) = -ilm{G;;(R, R, w)} - the local den- 
sity of states (LDOS). In particular, in the absence of 
the adatom, the LDOS becomes the unperturbed surface 



2. Calculating pa{R.,i^) for the subsystem adatom+surface 

To determine the LDOS in the presence of the 
adatom, we calculate GJ^(r3, rg, w) by Fourier transform- 
ing G;,3k^_,(w), 

G;(r3,r^,^) = e'^---^e-^^'r'^'^Gl^^,Juj), (26) 

kak;, 

assuming there is no tip-to-surface coupling, i.e., consid- 
ering H = H2 + H3 + H23- Following the procedure we 
described in the previous section - equation of motion + 
analytical continuation - we obtain^^ 

k2k2 

5k3.(^).9k^a(^)GL.k^,.(^), (27) 

where 51,^^(0;) = {uj^e^.+iS)-^ and ^^^k^ ^l^w^ai"^) = 
G22o.(w) is the adatom retarded Green function. Using 
Eq. (27) in Eq. (26) we find 

ikalra-rg) 

G;(r3,r^,c.) = Y^——-^ + |i23p X 



k; 

jka-rg 



^ w - efe3 + Z(5 ^ a; - e^/ -f- id 



G^2.(^)-(28) 



Let 



E — 



= i^(r3,w)+^/(r3,w), (29) 



where R{r3, ut) and /(r3, w) denote the corresponding real 
and imaginary parts given by 



R{r3,io) = pt / dx- Jo (kFr^VTTx) 

J-i {%-x) +6^ 



and 



(30) 



/(r3,c.) = -7:ppo (^kpr^^l + ^ , (31) 



5 



with kp being the Fermi wave number. We ean then 
write the LDOS as 



in this equation via Eqs. (l)-(3) we find 



p,{R,Lo)^pUi + ^-^jnkFRji 



D 



[(l-<z2)ImG^2,-2gReGS2,M]|, 

(32) 

whereTg = 27r|t23|V3 and g = R{R,uj)/ I{R,uj). In order 
to determine the transmission coefheient, we only have to 
calculate the adatom retarded Green function 6220- (w), 
obtained here using the Hubbard I approximation.'^^ This 
approximation accounts for the Coulomb interaction and 
consists in factorizing the higher-order correlation func- 
tions appearing in the resulting equation of motion for 
^22(t('^)- a result we have'^^ 



1 



with 



5L(^) = 



uj - e^y -U{1 - fig) 
{uj - e^){uj - e„ -U)' 



(33) 



(34) 



where a = — cr, rig. is the average occupation and Y7 is the 
self energy related to the coupling between the adatom 
and the host surface, T,^ = — -IFa. 



B. Resonant + non-resonant transport 

In this section we determine the total current - tip-to- 
adatom-to-host -|- tip-to-host - flowing in the system. In 
addition to the non-resonant (tip-to-host) current, cal- 
culated in Section A, here we consider the contribution 
from the direct tunneling of electrons between the tip 
and the adatom (resonant current). The corresponding 
additional term H12 is now taken into account and the 
Hamiltonian describing the system is given by Eq. (1). 
Note that for large enough R distances {R — >■ 00) we 
regain the non-resonant case. Here for convenience we 
perform the calculation in k space. 

We model the tip-surface, tip-adatom and adatom- 
surface couplings, respectively, by considering 



tl2 
tl3 



,0 --Ryi?,o 
''12'' I 



*23 = ^23; 



(35) 
(36) 
(37) 



where 1^2, ^13 and i23 are (constant) phenomenological 
parameters and Rq gives the exponential spatial decay 
for the coupling between the tip and the adatom as the 
tip moves away from it. 

The current flowing into the tip or leaving it can be 
determined from Eq. (5). Calculating the commutator 



3 

= (-^lj4iaCk,a + iyc[^.^Ckia) (38) 



3=2 ki.kj 

Substituting this result into Eq. (5) we obtain 

3 

= -«eZ Z {-*ij(4iaCk,T) +tlj(4^^ck,^)| 



j=2 ki,kj 



(39) 



or in terms of the lesser Green function G^.|^^ cr(^ii)i 



/r = 2eRe <^ ^ ^ h,G<^^^,^{t,t) \ , (40) 

[j=2ki,k, J 

where G^,]^^^^{t,t) = i{cl.^„{t)c\^^„{t)) . Equation (40) is 

equivalent to Eq. (11) when ti2 = 0. Now we must find 
an expression for the lesser Green function. Following the 
procedure described in section A, below we write down 
the equation of motion for the contour-ordered Green 
hmction G]^^\^^^„{t,t') = -i(T'cCkj.a(T)4^^(T')) 



d 



dr 



7 +eki ) Gk,ki,ff (r,r') = - ^^tliGi,^-k^.AT,T'), 

(41) 



;=2 k, 



or in the integral form 

3 



Gk.ki, ^(t,t') = ZZ*i' / ^^i'^k,k,,<T(T, ri)gki<T(n,T'), 



1=2 k 



(42) 

where gkicr('''i, t') is the tip free-electron Green func- 
tion, and then analytically continue Eq. (42) to find 
^kjki a(*'^')- Using this result in Eq. (40) we have 



II = 2eRe j ^ ^ iyi^, / dt^ 
U,;=2ki,k,,k, 



(43) 



Taking the Fourier transform of the above expression we 
find 

T, '"'i2K,k;,,(«^)9f,,(«^) + G£.ki,,('^)9i;„(")l + 

kik2k2 

E ^I2it3k3 [Gi;.k3..(^)5^,.(^) + Gk,k3,.(^)<ff£,.(c^)] - 

kik2k3 
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kikako 



kikak;, J 

(44) 

Using Eqs. (35)- (36) we can rewrite Eq. (44) as 

\t%?[Ghagt + G%,gl,]^, (45) 
where we have introduced the definitions 



= E s'^^^^Gkakj^cr, 
k3kj 

^iScr = E 6~''^^^Gkjk3,a-i 
kjk3 

G33.= E^'^'^"'^^''Gk3ki. 



(49) 
(50) 
(51) 



9la (^) 


= E^ki^^^)' 

ki 


'-^32(7 


= ^ e^i^-^G 

k3k2 


'-^23cr 


k2k3 


'^33a 


^ E 6*'''"'''=' 
k3k3 



2k3,0-' 



(46) 

Note that from Eq. (44) we regain Eq. (21) in the hmit 
^12 = 0, i.e., when the tip is far away from the adatom. 



1. Matrix Green function formulation 

We can see from Eq. (44) that G^j^.^ ^{t, t') is coupled 
to other Green functions. In order to find these Green 
functions, we have to apply the equation of motion tech- 
nique to the corresponding contour-ordered Green func- 
tion for each one of them and then perform an analyti- 
cal continuation to obtain the respective G^'''. After a 
straightforward calculation we find 



Gkikj,<T(T, r') = (5k,kj,9kjo-(T,T') + 

E / f^^Gk„k,,.(T,?)t*,.gk,<T(?,T'), (47) 



E 

l{l^3) kz 



i.e., a system of coupled equations for the Green func- 
tions. Here giija{T, t') is the free-electron Green function 
of the tip (j = 1), the adatom (j = 2) or the host sur- 
face (j — 3). These three Green functions g-kjcr{T,T') 
(j = 1,2,3) can be easily evaluated. Interestingly, by 
defining 



Gija — ^ ^ Gkikj,o-i 



(48) 



k3k!. 



where the sum is not taken over the spin indices, we can 
write down a Dyson equation of the form 

G.(t,t') =g.(T,T') + j driG,(T,Ti)Sg,(ri,T'), (52) 

with G„{t,t') being a matrix Green function whose ele- 
ments are defined following Eqs. (48)-(51), i.e.. 



Gila 


Gi2ct 


Gl3cr 


G 


G2ICT 


G22CT 


G23cr 


G 


Gsia 


G32a 


G33cr 


G 


Gsia- 


G^2a 


G33cr 


G 



13o 



(53) 



Additionally the self-energy is given by 



/ ti2 t% \ 

t*i2 t23 

^23 

V*;; y 



(54) 



and 



Sa 



/ 9ia 

92a 

gsa 93- a 

V 53^ 93a 



with the matrix elements'*^ 

9j<yiT,T') = 2^ ,9k,a(T,T'), j = 1,2,3 
k, 

53.(r,r') = Ee*''"V3.(^,^'), 

k3 

ff3.-(^.^') = E^"''"''5k3.(T,r'). 



(55) 

(56) 
(57) 
(58) 



Performing an analytic continuation in Eq. (52) we 
obtain the Dyson equation for the retarded Green func- 
tion 



and the Kcldysh'*'^ equation 

G< = G;g;"g<g^"G^, 

where 



9io 







r < 

= 



9'2a 



\ 



r,< r,< 



kik. 



53." 53 
V g;'< f,'< ) 



(59) 
(60) 

(61) 
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The advanced Green function g° is given by g° = [g^]*. 
From Eqs. (59) and (60) we see that if g^ and g^ are 
known wc can determine immediately and G^, and 
so the spin-resolved current, Eq. (45). The first nonzero 
(diagonal) elements {gf^ and gl^) in Eq. 61 are 



9la 



{lo) 



ki 



(i^) = ^ 2Trifi{uj)S{uj - ekio-) 

ki 



ki 



P 



Ekicr 



(62) 



inS{uj - ekiff) 



C. Parameters and a summary of the numerical 
technique 

Our main task is to determine the current from Eq. 
(45). To this end, we first calculate GJ^ and G^ from 
Eqs. (59) and (60), respectively. Then we substitute the 
relevant matrix elements in Eq. (45). Note that in the 
presence of the Coulomb interaction, g2^{ui) and g2„{uj) 
depend on the adatom occupation rig.; so do GJJ. and 
G^ . This implies a self-consistent calculation, where no- 
is calculated iteratively via 



= Ai(a;) - i7r/9icr(w), (63) 
where fi{uj) is the tip Fermi distribution function, pia{i^) 
is the tip density of states and Ai(a;) = -PX^ki (zT^ei; — ) ' 
where P stands for the Cauchy Principal Value. The re- 
tarded adatom Green function .920- ('^) given by Eq. 
(34). The lesser component can be calculated straight- 



forwardly from the relation g2^{uj) - 
Acr{uj) — 27r(l — ng-)(5(aj — e^) + 2Trng-S{u! - 



jA(j(w), where 
-[/), (64) 



and is the average spin-resolved occupation of the 
adatom. The third diagonal element of g*" is given by 

1 



k3 k3 

= A3(w) - j7rp3(a;). 



P 



uj - 



— ITT (5 (a; — Ekg) 
(65) 



where ^3(0;) = is the 2D density of states of the 
surface defined at the end of Sec. II A 1, and A3 (a;) = 



P 



Ska (l 



For the corresponding lesser Green 



function wc find 



9^0 



E.9k3.H 

ka 



2^«/3(^)P§, 



(66) 



where f 3(1-0) is the Fermi distribution function of the host 
surface. 

Finally, we should calculate the off-diagonal elements 
of the matrix g^. For the retarded Green function we 
have 



.93, 



.ika-R 



dw 
27ri 



G2<2o(^). 



(69) 



As a matter of simplification we use the wide-band 
limit for the tip, so the density of states picr(w) is 
taken as constant pi (pg is already a constant), eval- 
uated at the Fermi level. This is a good approxima- 
tion when eV, /cbT << D, where D is the band half 
width. The ferromagnetism of the tip is introduced via 
the density of states pia = Pi{l ± p), where p is the 
tip polarization and the -I- and — signs apply to spin up 
and down, respectively.'^^ Since the characteristic tun- 
neling rate between the tip and the adatom is given by 
Tia - 27r|ii2|2pio, wc find Fi^ = 2n\h2\^p1{l±p), which 
is the standard phenomenology to account for the fer- 
romagnetism of the electrode.'*'* Analogously, the tun- 
neling rate between the adatom and the host surface 
is Fg = 27r|i23pP3 = ^o- In our calculation we take 
Fo = 10/ieV as the energy scale. All the phenomeno- 
logical parameters used in this paper are summarized in 
tabic I. 



Parameter 



Magnitude 



Band half width 
Adatom Charging Energy 
Tip chemical potential 
Host chemical potential 
System temperature 
Tip degree of polarization 
Decaying factor of ti2 



D = lOOOFo 
U = SOPo 
^1 = -eV/2 
P3 = eV/2 
ksT = To 
p = 0.4 
-Ro = l/kp 



TABLE I. Parameters used in the self-consistent calculation. 



1 yo 



D 



X)2 + 



Jo{kFRvT+x) 



inp^Jo 

For the lesser Green function we find 



93o 



2^ip^/3(cj)Jo kFRjl + 



kFRjl + -]. (67) 



(68) 



The Green functions g^^^ and g^^^ have exactly the 
same expressions as Eqs. (67) and (68), respectively. 



Note that in Table I we define the origin of the en- 
ergy scale such that /ii = //3 = eo- = at zero bias, i.e., 
the adatom energy level is aligned with the chemical po- 
tential of the leads in the absence of an applied bias. In 
order to be consistent with Fq ~ 27r 1^23^/33 and the value 
of pg for typical semiconductors, e.g, GaAs, the factor <23 
is taken as To/VTO. We assume t^r^ = O.OK23 and adopt 
values for and pi consistent with Fq = 2n\ti2\^Pi- 
Note that for eV = ±30Fo the adatom can be occupied 
by a single electron since is within the conduction win- 
dow (the energy range between pi and p^) and + U is 
without this range. On the other hand for eV = ±150Fo 
the dot can be doubly occupied since both £„■ and + U 
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lie inside the conduction window. The parameter Rq con- 
trols how fast the coupling ti2 decays in space when the 
tip moves away from the adatom. We take it equal to 
kp^. Hereafter kp^ will be used length scale. 



III. RESULTS 

A. Single Occupancy 

Figure 2(a) shows the adatom occupations against the 
lateral distance R between the tip and the adatom. We 
study both forward {eV = +30ro) and reverse {eV = 
— SOFo) bias cases. As mentioned before, for \eV\ = SOFq 
the adatom cannot be doubly occupied, since + U lies 
above the Fermi energy of the source. For V > the host 
NM surface is the source and the FM tip is the drain of 
electrons, i.e., the electrons flow from the NM surface (+ 
adatom) to the tip. For < we have the opposite. 



3 

o 
o 

O 



o 

-a 

< 



0.6 
0.4 
0.2 


-0.2 
-0.4 



(a) 




thicker solid lines: spin up 




thinner solid lines: spin down 




solid circles: magnetization 


• 

— • 

1 1 


black curves: V>0 
grey curves: V<0 
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0.0067 



0.0066 
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FIG. 2. (a) Adatom occupations and magnetization and (b) 
spin-resolved currents against the tip position R. Both neg- 
ative (grey line) and positive (black line) bias voltages are 
shown. For \eV\ — SOFo only single occupancy (i.e., e„ + U is 
much higher than the emitter chemical potential) is allowed. 
This results in the spin-diode elfect: for negative bias volt- 
ages the current is polarized for all R values while for positive 
biases the current polarization goes to zero when the tip is 
close to the adatom. Inset: zoom of each spin component 
of the current in the range where the tip-adatom coupling is 
negligible. In (c) we also show the LDOS (thick solid black 
line) evaluated a,t ep. The LDOS displays Friedel oscillations 
which show up in the polarized current. The vertical scale for 
the LDOS is not shown. 



For R = and > (black lines) we find < nj,. 
This is reasonable since spin up electrons tunnel from the 
adatom into the tip much more easily than the spin down 
ones, due to the larger spin up density of states, pia = 
Pi{l ± p) {p > 0), which implies Fi-f- > F14,, and gives 
rise to a larger spin down population. For a constant 
bias voltage, as the tip moves away from the adatom [see 
Eq. (35)] the tunneling rates Fi-j- and Ti^ decrease, but 
the incoming rates F3 stay the same, this results in an 
increase of both rif and n^. We note, however, that the 
spin down population (thinner solid black line) increases 
more quickly than the spin up one (thicker solid black 
line), this is because Fi^ < Fi^, see Fig. 2(a). As the rif 
adatom population increases with R, the ni one tends 
to be more blocked due to the spin-dependent Coulomb 
blockade. The interplay between the Coulomb blockade 
and the decrease of the tunneling rates Fi-|-^^ makes 
reach a maximum, subsequently decreasing to attain the 
limit n-f ~ ni ~ 0.5 for large enough i?'s. 

In contrast, for eV = — BOFq, n-^ > for small R val- 
ues, see in Fig. 2(a) the solid thicker and thinner grey 
lines. This is a consequence of Fi^ > Fi^, which means 
that more spins up tunnel to the adatom. Besides, the 
outgoing rates F3 (the same for the up and down compo- 
nents) is smaller than Fi-j-, which results in a larger spin 
up accumulation in the adatom. As the tip moves away 
from the adatom, Fi-|- and Fij^ go to zero exponentially 
and the populations n-f- and are completely drained out 
into the host surface, thus resulting in an empty adatom. 

The magnetization m = — is also shown in Fig. 
2(a) (solid circles). Observe that for small R the adatom 
is spin down polarized for V > and spin up polarized 
for V < 0. As R increases, m tends to zero for both 
positive and negative bias voltages. However, m tends 
to zero much slower for V > than for y < 0, a conse- 
quence of the interplay between the Coulomb interaction 
(spin-Coulomb blockade) and the tunneling rates Fio-, 
that change with the tip position as it moves away from 
the adatom. 

In Fig. 2(b) we present the spin-resolved currents for 
both eV — ±30Fo. The spin-diode effcct^'^ can be clearly 
seen for small values of R. While for y > (black lines) 
we find w /j. for small R values, ioi V < (grey lines) 
we observe > Ii- This shows that the current polariza- 
tion can be controlled via both the bias sign and the tip 
position. In the case oiV<0, we have Fi-|- > Fa > Fij,, 
i.e., the spin up population is greater than the spin down 
one, m > 0. As a consequence, in the absence of the 
Coulomb interaction in the adatom, > (the case 
U = resembles the curves in the double occupancy 
regime {eV >> U), see Fig. 3. In the presence of U, 

is suppressed, since tends to be more blocked than 
[see Fig. 2(a)], which results in an enhancement in 
the difference between If and I^. For I^ > 0, the mag- 
netization changes sign m < 0, now the spin up popula- 
tion tends to be more blocked, and If is more strongly 
suppressed compared to If, interestingly attaining values 
close to I^. The amplification of If compared to If for 
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V < 0, when the tip is closer to the adatom, does not 
occur in the double occupancy regime {eV = ztlSOFo) as 
we will see in the next section. 

In Fig. 2(c)-(f) we show the current for a range of R 
in which only the direct tip-host tunneling (non-resonant 
transport) is relevant. Note that and /j, tend to dis- 
tinct plateaus for large enough i?'s. These plateaus corre- 
spond to the background current between the tip and the 
host surface, given approximately by Eq. (25). By com- 
parison with Eq. (45) we plot in dashed line the current 
obtained via Eq. (23). In the large- i? limit, we expect an 
agreement between both equations, since Eq. (23) was 
derived in the case of negligible tunneling between the 
tip and the adatom (see the solid black and grey lines). 
The minor difference between the two results is due to 
Eq. (23) having been obtained in the limit of small tip- 
surface coupling parameter t^^-^^ The LDOS evaluated 
at the Fermi level, pa{R,0), is also shown in Fig. 2(c); 
it oscillates around the unperturbed surface density of 
states pg. Friedel-like oscillations are seen for both spin 
components, thus reflecting the oscillations in the LDOS 
due to the scattering center (adatom). Note that Friedel 
oscillations have been seen experimentally in a variety of 
systems. 
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B. Double Occupancy 



FIG. 3. Similar to Fig. 2 except that \eV\ = ISOFo here. This 
bias allows double occupancy of the adatom. In this regime no 
spin-diode effect is observed. The spin polarized currents are 
the same (in modulus) for both positive and negative eV. For 
forward voltages (black curves) the adatom becomes doubly 
occupied for large R values (n-|- + n_i = 2). 



Figure 3 shows the spin-resolved (a) adatom occupa- 
tions and (b) currents in the double occupancy regime, 
i.e., when the bias voltage is large enough {\eV\ = 150Fo) 
to allow for two electrons of opposite spins in the adatom 
at the same time. For V > (black lines), as R increases 
both and ni populations tend to unity and the adatom 
becomes doubly occupied (n^ + = 2). This is so be- 
cause electrons can jump into the adatom but cannot 
leave it for large R values. In contrast, for V < both 
and vanish as R increases because the electron 
source (tip) decouples from the adatom. Observe also 
that the magnetization m is enhanced as R tends to zero 
and switches sign depending on whether eV is positive 
or negative. 

The current in the double occupancy regime [Fig. 3(b)] 
has a similar behavior for both positive and negative bi- 
ases. Note that |/-t-| > |/4,| for F > and ^ < 0, in 
contrast to the single occupancy regime where we find 

~ for > [Fig. 2(b)]; hence no spin-diode ef- 
fect is observed here. In Figs. 3(c)-(f) we show a zoom 
of the spin-resolved currents in the range of negligible 
tip-adatom coupling. As before, we observe Friedel oscil- 
lations, which rcfiect the disturbance in the LDOS due to 
the localized impurity (adatom). The dashed black lines 
in Figs. 3(c)-(f) show the current obtained via Eq. (23). 



C. Current Polarization 



Figure 4 displays the current polarization, p = {Ij — 
+ as a function of R. Both single {\eV\ = 
30ro) and double {\eV\ = 150ro) occupancies are shown. 
For eV — — 30ro (solid grey line) the polarization is ap- 
proximately constant ^ 40% while for eV = -t-30ro (solid 
black line) it is strongly suppressed when the tip is close 
to the adatom. This is a type of spatially resolved spin- 
diode^"^ , which allows a polarized current to flow only for 
reverse bias. In the case of double occupancy, though, 
both positive and negative biases present a similar be- 
havior with a 40% current polarization away from the 
adatom and a slight suppression as the tip moves closer 
to it. This 40% current polarization for large i?'s in 
the single- and double-occupancy cases follows straight- 
forwardly by calculating p using the spin-resolved non- 
resonant currents in Eq. (25). The inset shows a blow 
up of the current polarization and also reveals Friedel 
oscillations. 
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FIG. 4. Current polarization against R. In the single occu- 
pancy regime (eV = iSOFo) the current polarization is ap- 
proximately constant for reverse bias {V < 0) and goes down 
to zero for direct bias {V > 0) when the tip is close to the 
adatom. In the case of double occupancy {eV — ilSOFo) the 
current polarization is suppressed for small R, however it does 
not vanish; notice that the dotted lines lie essentially on top 
of each other. For large _R's, the polarization for both single 
and double occupancies tends to a plateau and presents tiny 
oscillations (inset). 



surface with a single adsorbed atom. Due to Coulomb 
interaction in the adatom the system can operate as a 
spin-diode when the tip is nearby the adatom. In the 
singly occupied case and direct bias (T^ > 0), the current 
polarization can vary from zero up to 40% depending on 
the tip position. For reverse bias, though, the polariza- 
tion is pinned close to 40% for all tip positions. In the 
double occupancy regime the current polarization is the 
same for both forward and reverse biases, with a slight 
suppression as the tip moves closer to the adatom. Ad- 
ditionally, the adatom magnetization can be tuned by 
varying the tip position and its sign can switch depend- 
ing on the bias. Finally, we have also found spin-resolved 
Friedel oscillations in the current as the tip moves later- 
ally away from the adatom, thus reflecting the oscillations 
in the surface LDOS induced by the adatom acting like 
an effective impurity. 
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